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Abstract 

We present calculations of the tunneling density of states in an anisotropically paired superconductor 
for two different sample geometries: a semi-infinite system with a single specular wall, and a slab of 
finite thickness and infinite lateral extent. In both cases we are interested in the effects of surface 
pair breaking on the tunneling spectrum. We take the stable bulk phase to be of d^2_y2 symmetry. 
Our calculations are performed within two different band structure environments: an isotropic 
cylindrical Fermi surface with a bulk order parameter of the form A ~ fc^ — fc^, and a nontrivial 
tight-binding Fermi surface with the order parameter structure coming from an anti-ferromagnetic 
spin-fiuctuation model. In each case we find additional structures in the energy spectrum coming 
from the surface layer. These structures are sensitive to the orientation of the surface with respect 
to the crystal lattice, and have their origins in the detailed form of the momentum and spatial 
dependence of the order parameter. By means of tunneling spectroscopy, one can obtain information 
on both the anisotropy of the energy gap, lA(p)l, as well as on the phase of the order parameter, 
A(p) = |A(p)|e'*'(P). 
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I. INTRODUCTION 

In this paper we study the density of states spectrum of a d-wave superconductor in the vicinity of a specularly re- 
flecting surface. The anomalous spectral properties of anisotropic superconductors and superfluids have been discussed 
in several publications p]-p^. The present work lo gica lly follows our examination of the system's order parameter 
and surface free energy which was presented in Ref. |13| (from here on we refer to this paper simply as [I]). Together, 
these studies are part of an ongoing effort to systematically examine measurable features which are typical of d-wave 
superconductors. We hope, as before, to contribute to the fund of theory necessary for their definitive experimental 
identification. 

As discussed in [I], one expects that an anisotropic order parameter will respond strongly to pair-breaking effects 
near a wall which is, again, why we have chosen this setting as our arena of study. Moreover, the resulting spectrum of 
excitation states is quite sensitive to the details of the entire system and is also accessible to a variety of experimental 
probes. Accordingly, there are any number of spectral dependencies one could examine which may well yield significant 
clues as to the true symmetry of the order parameter. We have chosen three for this paper: 1) surface to lattice 
orientation, 2) realistic Fermi-surface geometries, and 3) finite size effects. The results presented here are not in any 
way intended as exhaustive studies of any of these areas. Instead, we selected out of a rather large collection of 
calculations a condensed subset which we felt was either representative of general behavior or exceptionally suggestive 
of potential experimental verification. 

We employ the quasiclassical formulation of superconductivity for the calculations of this paper and refer the 
reader to [I] for a more detailed discussion of these methods. The equations of particular importance for the current 
work are collected in section II, and the results are displayed in the two subsections of section III. In section IV we 
discuss various aspects of the physical significance of our resulting data and we conclude with a summary of the most 
prominent points. 
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II. QUASICLASSICAL THEORY 



A. Basic Equations 



In [I] we presented the fundamental equations of the quasiclassical theory in the imaginary-energy (Matsubara) 
representation. These considerations provide a sufficient framework for the study of the thermodynamic properties 
of the system. A discussion of the spectral properties of the system, however, requires the real-energy version of the 
Green function equations presented in [I]. We therefore must compute the retarded propagator, g^, which can be 
defined in terms of the Matsubara propagator as follows (see, for example, Ref. [Q): 

(p/, R; e) = 5 (p/, R; e„) |,,,,_,+,^ . (1) 
The 2x2 matrix propagator g^ obeys a transport equation of the same form as in the case of imaginary energies: 

ef3-A(p/,R) ,g«(p/,R;e)] iv/ (p/) • Vr5^- (p/, R; e) = , (2) 
with the same normalization condition, 

[5^(p;,R;6)]' = -7rM. (3) 

Here -v f (pf) is the momentum dependent Fermi velocity, and A (p/, R) is the superconducting order parameter (or 
pairing self-energy). For anisotropic superconductors A(py,R) exhibits, in general, a nontrivial dependence on the 
momentum, p^, and the position, R, which must be calculated self-consistently from the gap-equation. The method 
for carrying out the self-consistent solution for A was discussed in detail in [I] and will not be repeated here. Instead, 
we simply carry over the results obtained in [I] for the self-consistent order parameter and use them as input in the 
solution of equations (||) and (^ for g^ . 

The propagator g^ carries detailed information on the quasiparticle excitations. In the quasiclassical scheme, a 
quasiparticle state is characterized by its excitation energy, e, measured from the Fermi energy, and by a Fermi surface 
coordinate, p/, which is the projection of the quasiparticle momentum onto the Fermi surface. The total number of 
quasiparticle states (per spin) with energy, e, in the energy interval c?e, and Fermi momentum, p/, in the interval dbpf 
is given by 

dA^(p/,R;e) = iV/n(p/)i.(p/,R;e)dp/de, (4) 

where Nf is the total density of states in the normal state as measured, for example, by specific heat experiments. 
Band-structure effects lead to the anisotropy factor 7i(p/), which is normalized such that § dpf n{pf) = 1. Finally, 
the effects of superconductivity on the density of states are carried by the dimensionless quantity i^(p/, R; e), 

(p/, R; e) = {Tr [r, g^ (p^, R; e)] } , (5) 

which we refer to in the remainder of the text as the "superconducting DOS factor" . In the normal state, i^(p/, R; e) — 
1. In what follows we focus on the calculation of the quantity i^(p/, R; e), and, in particular, how its qualitative features 
are affected by the presence of a wall. 

We consider two different sample geometries: a semi-infinite system occupying the space x > 0, and a slab of finite 
thickness, w, and infinite lateral extent oriented perpendicular to the a;-axis. In both cases we consider perfectly 
reflecting (specular) walls. Continuity of the propagator along a classical trajectory requires | p5| , p^ : 

■9^(P/,in'I^wall!^) = 3^(P/.out' R-waU^ (6) 

which yields the appropriate quasiclassical boundary condition for a specular wall. An in-coming trajectory is defined 
as one in which vj(pj ■^^) ■ n < while an out-going trajectory requires that vj(p^ out) ■ > (here h is taken to be 
the surface normal pointing into the sample). The out-going momentum vector, Pjout' then given in terms of the 
in-coming momentum vector, p^ jj^, by requiring that the momentum component parallel to the surface be conserved 
upon reflection, 

Pj,out = Pjin- (7) 
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Note that energy conservation is imphed in equation (|[) since e^^ = Cq^^. 

We restrict our attention to the case of simple Fermi surfaces in which any given "in-coming" trajectory is uniquely 
connected to an "out-going" trajectory via equation (|^). This boundary condition would need to be modified for 
systems with more complicated Fermi surfaces (e.g. multiple Fermi sheets, etc.), and additional phenomenological 
parameters would be required to determine the probability of a given in-coming trajectory being reflected into any 
one of the multiple out-going channels. 



B. Nontrivial Band-structure 

Before the self-consistent order parameter can be calculated and a complete solution of equation (||) can be carried 
out, we must first specify the form of the Fermi surface, the corresponding Fermi velocity, v/(p/), the anisotropy 

factor, n(p/), and the quasiparticle pairing interaction, V (pf, p^-^ (see equation (4) in [I]). We approach this problem 

in two different ways. First we consider an isotropic cylindrical Fermi surface with a pairing interaction that can be 
parameterized in terms of simple trigonometric functions (as in [I]). The second approach, which we carry over from 
Refs. and |8|, is to use actual tight-binding fits to the band-structure of high- Tc materials, along with a pairing 
interaction derived from a microscopic spin-fiuctuation theory. The first of these two approaches was discussed in 
detail in [I], and we therefore focus on the second approach in what follows. 

Following Radtke and Norman , we approximate the band-structure of the high-Tc cuprate superconductors by 
a 2d tight-binding model of the form: 

E{p) = to + — [cOSpx + cos Py] + t2 cos Pa; COSpj, -|- — [cOS 2px + COS 2py] 

+ ^ [cos 2px cos Py + COS 2py cos Pa,] + cos 2p,j. cos 2py. (8) 

Given a Fermi energy, Ef, we use the relations Ef — E{pf) and = WpE{p)\Ef to extract the Fermi wavevector, 
P/, at every point on the Fermi surface, along with the corresponding Fermi velocity, v^(p/). These quantities then 
allow us to compute the anisotropy factor, n{pf), in the usual way. Once the band structure data is known, the 
pairing interaction kernel can be evaluated using the weak-coupling version of the spin-fiuctuation model of Radtke 
et. al. as follows: 

^(P/,P/) =-^o^(P/-P/), (9) 

where 



(10) 



1 + Jo (cos (7a, a -I- cos qy a) 

One obtains this weak-coupling version of the interaction presented in Ref. \V7\ by integrating-out the high-energy 



degrees of freedom |14 . In doing so one introduces a cutoff, £co, and a renormalized interaction parameter, Vq, 
which differs, in general, from the bare interaction parameter of Ref. [ pT| . The renormalized interaction strength, 
together with the energy cut-off, can always be eliminated in favor of the superconducting transition temperature, 
Tc- Hence our pairing interaction is determined by two phenomenological parameters: the transition temperature, 
Tc, and the exchange coupling parameter, Jq. The transition temperature is written in terms of the energy cut-off (in 
the weak-coupling theory) via: 

T,^ = l.lSecoe-^/^" (11) 

where the are dimensionless coupling constants in the various symmetry channels (see below). 

Once the quasiparticle pairing interaction is known, the superconducting order parameter can be computed by 
self-consistently solving the weak-coupling gap equation (equation (4) in [I]) in the Matsubara representation. This 
equation can be recast into a more convenient form by introducing a complete set of eigenvalues and eigenfunctions 
of the interaction kernel: 

A^ry^ (P/) = f dp'f nip'f) V (p;, p}) r?^ (p'^) , (12) 
where the eigenfunctions, 77^ (p/), satisfy the orthonormality condition: 
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^dp/ ?i(p/)?7^ (p/)77* (p/) = (5^^. (13) 

If we now express the superconducting order parameter as a linear combination of these eigenfunctions, 

A(p/,R) = ^A^(R),7^(p;), (14) 

then the gap equation may be rewritten as a series of equations for the individual gap amplitudes as follows: 

(R) = X^Tf^<f dpf n{pf) rj; (pf) f {pf,R; e„) . (15) 

The eigenfunction, 77^ (p/), associated with the largest (attractive) eigenvalue determines the symmetry of the bulk 
stable phase, at least near Tc. We assume that this symmetry persists to all temperatures. The effects of the other 
symmetry components however (whether they are associated with attractive or repulsive eigenvalues), may become 
noticeable in the vicinity of surfaces ||l3], defects, in the cores of vortices, etc. 





Fig. 1. Band-structure data: (a) The Fermi surface in the first Brillouin zone. We parameterize this surface by 
an angle (j> = [0, 2-k] taking the point indicated by the solid dot as the origin. In the repeated zone scheme our 
Fermi surface is then a single continuous line, around the origin, which is single- valued in </!>. (b) The a;-component 
of the Fermi velocity normalized so that the maximum is 1. (c) The anisotropy factor of the normal state density 
of states. The dotted lines show the corresponding values for the isotropic cylindrical Fermi surface. 



In practice, equation ( [l5| ) must be solved numerically, which first requires the numerical evaluation of the eigen- 
values, A^, and and eigenfunctions, 77^ (p/), as defined in equation (^2|). Since we would like to expand our order 
parameter in terms of eigenfunctions of well defined symmetry, we first project out the pieces of the interaction 

matrix, V ^p/,Py^, which transform like each of the four one dimensional representations of the D^h (tetragonal) 
group (i.e. Vai ^p/,p'^^, (pfiP'f^i (pfiP'f^i (pfiP'f^)- We then use each of these matrices in place 

of the V ^p/,p'^^ in equation (|l^ ) to solve for the eigenvalues and eigenfunctions corresponding to each symmetry 

classification. These eigenfunctions then allow us to construct a gap consisting of specific symmetry components, 
while the eigenvalues give the coupling strengths in the various symmetry channels. 



We calculate the Fermi surface data using the parameter values reported in Ref. [ 
(i.e. to = 0.0, h = -1.0, t2 = 0.38, = -0.09, t4 = 0.0, <5 = 0.0, aU in units of eV). T 



_J for the YBC0(2,bf) sample 
le Van Hove singularity for this 
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system lies at about E = —AGeV |19), and we take a Fermi energy slightly above this value at Ef = —.30eV. The 
resulting form of the Fermi surface in momentum space is shown in Fig. la. Repeating the Brillouin zone (repeated 
zone scheme) allows us to parameterize our Fermi surface by a polar angle (j> around the k=(— 7r/2,— 7r/2) point in 
the Brillouin zone (denoted by a solid black dot in Fig. la). In Figs, lb and Ic we plot the Fermi velocity, v/(0), and 
the normalized density of states anisotropy factor, n(0), as functions of (p. The factor n((j)) is obtained from equation 
(4), yielding 

(27r)^ |p/(0) • V/((/))| 

and is normalized such that d(j)n((j)) = 1. 




Fig. 2. The momentum space dependence of the Bi and A2 order parameter components as computed from the 
weak-coupling spin-fluctuation model discussed in section II-B. The dotted hues show the corresponding results 
when a cylindrical Fermi surface is used. 



Once the band-structure data is known, we may then proceed to compute the spin-fluctuation pairing interaction 
via equation (|^). Choosing the dimensionless parameters Vq = —1.0 and Jq = 0.4, we find two attractive symmetry 
channels: Bi and A2. The Bi channel has the highest transition temperature while the A2 channel has a Tc that 
is reduced by about 3/7 (the transition temperatures for the repulsive symmetry channels are set to zero in what 
follows). Within the dominant symmetry channel, the next to leading coupling constant is more than two orders 
of magnitude smaller than the leading coupling constant and therefore only the leading contribution is considered. 
The angular distribution of the order parameter components corresponding to the Bi and A2 symmetry channels are 
depicted in Figs. 2a and 2b. 



III. RESULTS 



We consider in both of the following subsections a surface geometry of the same form as in [I] . The normal to the 
surface is taken to lie along the ii-axis, which need not, in general, be parallel to the crystal lattice basis vector a, 
from which we measure all our angles. We denote the angle between a and the surface normal, x, by (jjo, and the angle 
between a and the momentum of an out-going trajectory, Pj^out' ^- This geometry is represented graphically in 
Fig. 1 of [I]. We denote the transition temperature for the bulk symmetry channel {Bi) by Tc while the transition 
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temperatures for the other symmetry channels are denoted exphcitly by TcAi, TcA2i TcB2- All length scales are 
measured in units of the coherence length which we define in the usual way by ^ = f\max / ''^ ^max~buik- 

A. Semi-Infinite System 

We consider a semi-infinite system occupying the space a; > with a perfectly reflecting (specular) wall located 
at a; = 0. The order parameter and free energy of this system were calculated in detail in [I] (considering an 
ideal cylindrical Fermi surface) for several different symmetry parameterizations of the order parameter. From our 
calculations of these quantities in the case of the tight-binding band-structure and spin-fluctuation pairing interaction 
discussed in section II-B, we conclude that the thermodynamic properties of this system are quite insensitive to the 
details of the Fermi surface geometry, as long as the band-structure does not become too pathological. Hence, we 
do not discuss the stability of the self-consistent order parameter calculations here. Instead, we use these results to 
calculate g^' from equation (||), and then focus our attention on the excitation spectrum. This includes the response 
of the excitation spectrum to an admixture of subdominant components and to nontrivial band-structure effects. 





p>0^ 




p-0 




> 




X 

\ V, . 


vacuum 


SC 





Fig. 3. A typical quasiparticle trajectory reflecting from a specular wall in the tight-binding model discussed in 
section II-B. The in-coming leg of the trajectory is indicated by a negative value for the trajectory parameter, p, 
with the point of reflection being at p = 0. Note that for an anisotropic Fermi surface the out-going angle is not, 
in general, equal to the in-coming angle. 

In Fig. 4 we show some representative curves for the superconducting DOS factor (defined in equation (|^)) at x = 0, 
together with profiles of the order parameter along the corresponding "trajectories" (i.e. the trajectory of a classical 
particle approaching the surface with velocity v^(p^jj^) and leaving with velocity v/(pjQ^^)). These curves were 
calculated for the case of a cylindrical Fermi surface with an order parameter of the form A ~ fc^ — fc^ . A negative 
value for the trajectory parameter p corresponds to the in-coming segment of the trajectory, while the point p = 
corresponds to the point of reflection (a typical trajectory is depicted in Fig. 3). Note that a characteristic "edge" 
appears in the superconducting DOS factor corresponding to the magnitudes of the asymptotic values of the order 
parameter at p = ±cx). For the = Odeg and (/> = 45deg trajectories, these values coincide so that we find only a 
single feature in each case at e = [1 / 2) /S.max-buik and at e = [■\/i/2)/S.max-buik respectively. Both the (j) = 20deg and 
the <j) = 45 deg trajectories exhibit a single pole (or bound state) at e = 0. These states were discussed by Hu ||], and 
are a consequence of the Atiyah-Singer index theorem |^ . It should be emphasized that the presence of these states 
is a robust feature that will always occur (at e = 0) whenever the local order parameter (which is real in our gauge) 
changes sign asymptotically along a trajectory. The additional pole found just below the continuum in the = deg 
curve, on the other hand, can be interpreted as a bound state in the "potential well" formed by the order parameter 
along the = deg trajectory. This is not a robust feature, and will only occur if the well has sufficient depth and 
sufficient width. The location and weight of this feature therefore contain information on both the magnitude and 
the width of the order parameter suppression near the surface, which are, in turn, directly related to the degree of 
anisotropy in the gap. 
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Fig. 4. Order parameter and density of states in a semi- infinite system (cylindrical Fermi surface, Bi-symmetry): 
(a) the local order parameter for various (out-going) trajectory angles, (f>, and (b) the superconducting density 
of states for the same trajectory angles. The lattice to surface orientation is taken to be 0o = 30deg, and the 

temperature is T = OATc. 
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Fig. 5. Order parameter and density of states in a semi-infinite system (tight-binding Fermi surface, Bi -symmetry): 
(a) the local order parameter for various (out-going) trajectory angles, </>, and (b) the superconducting density of 
states for the same trajectory angles. The lattice to surface orientation is taken to be 0o = 30deg, and the 
temperature is T = 0.4Tc. 
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For essentially all surface to lattice orientations, the inclusion of nontrivial band-structure data results in rather 
minor effects on the density of states. In Fig. 5 we show the superconducting DOS factor for the same parameter 
set as in Fig. 4, only now we consider the band-structure data and pairing interaction discussed in section II-B. We 
have suppressed the presence of the attractive A2 symmetry channel by artificially setting its transition temperature 
to zero; thus Fig. 5 corresponds to a Bi-symmetry order parameter with a momentum dependence of the form given 
in Fig. 2a. Other than a redistribution of the location of the gap features, Fig. 5 is remarkably similar to Fig. 4. The 
fact that the gap features reorganize in this case is understandable from Fig. 2a where we see that the magnitude 
of the order parameter for most momentum angles has been slightly altered. However, the key point is that for a 
single-band Fermi surface of the form discussed in section II-B, we expect fairly little alteration in the results from 
the simple cylindrical Fermi surface model. Such Fermi surfaces have recently been obtained in connection with many 
of the high-Tc systems through the use of 2-dimcnsional tight-binding models and therefore may be considered 
as reasonable models for the band-structure in some of the cuprate systems. 
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Fig. 6. Order parameter and density of states in a semi-infinite system (cylindrical Fermi surface, B\ -\- B2- 
symmetry): (a) the local order parameter for various (out-going) trajectory angles, 4>, and (b) the superconducting 
density of states for the same trajectory angles. The lattice to surface orientation is taken to be cj>o — 30deg, and 
the temperature is T = 0.4rc. 



We have also made an extensive study on the effect of the mixing of sub-dominant symmetry components with 
the stable bulk (_Bi -symmetry) phase. We define a sub-dominant symmetry component as any attractive symmetry 
channel which has a lower transition temperature than the bulk symmetry channel. These effects were discussed in 
detail in [I] in connection with the calculation of the self-consistent order parameter and free energy. It was shown 
there that the most prominent effect comes from the mixing in of the B2 symmetry component since it allows for an 
effective rotation of the total order parameter towards its optimal (symmetric) configuration at the wall. This has a 
very pronounced qualitative effect on the structure of the order parameter near the wall, which is in turn discernible 
in the density of states. In Fig. 6 we plot the superconducting DOS factor, i/, for a system consisting of a Bi bulk 
phase coupled to a B2 symmetry component with a lower Tc. The choice of parameters is the same as in Fig. 4 for 
the purpose of comparison. The zero energy bound states are still present in the (f> = 20 deg and = 45 deg curves 
since the formation of these states depends only on the asymptotic values of A{p) which are not changed when surface 
mixing is present (i.e. the bulk phase is not changed). The pole in the </> = deg trajectory, however, has now vanished 
since the rotation of the total order parameter gives rise to a local maximum in the order parameter near the wall 
rather than a local minimum. On the other hand, an additional bound state now appears just below the continuum 
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in the (p = 45 deg curve. The appearance of this pole is a result of the broadening of the local minimum formed by 
the absolute value of the order parameter along the (/> = 45 deg trajectory when mixing is allowed to take place. 

So far our discussion of the tunneling density of states has focused on the value of the superconducting DOS factor, 
V, at the wall, since this is the case which is presumably accessible to tunneling experiments. For completeness, 
however, we show in Fig. 7 the behavior of v as we move into the sample. Here we have selected the <f> = Ab deg 
trajectory curve from Fig. 6 since it contains both bound state features. As we move away from the wall, v approaches 
its bulk form by continuously transferring weight from the bound states to the gap edge. The inset in Fig. 7 depicts 
the manner in which the bound state weights decay as we move into the bulk. While the zero-energy bound state is 
highly localized near to the wall, the upper bound state is spread out over a much broader region, which once again 
emphasizes the distinct nature of these two states. The ripples that develop in the superconducting DOS factor for 
e > ^max-buik ^rc Tomasch oscillations whose prominence is a consequence of the cleanliness of our system. 




0.0 0.5 1.0 1.5 



'max-bulk 



Fig. 7. The density of states for a _Bi-symmetry order parameter coupled to a sub-dominant i32-symmetry 
component as a function of distance from the wall. The inset shows the dependence on x of the weight of the two 
bound state delta functions (at e = 0.0 and at e = 0.83). 



The mixing in of an order parameter component transforming according to the A2 symmetry channel has also been 
considered. However, since the inclusion of an A2 component only allows for a distortion of the Bi order parameter 
in momentum space (rather than an overall rotation) , we find very little variation in the qualitative features of the 
density of states. 

We presented in this section our results for the angle-resolved density of states near a surface. We close this 
section by discussing the I-V characteristic for quasiparticle tunneling of an S-N junction where S refers to a d- 
wave superconductor. The differential tunneling conductance, dl/dV, is, at zero temperature, given by the following 



equation |23 



di 1 r^^ 

d<pD{P)iy{<t>,eV), (17) 

-Tr/2 



dV Rn 



where v{4),eV) is the normalized density of states (see equation (4)), and !?((/)) is the "barrier function" which is 
normalized such that 

Tx/2 

d(f,D{(f>) ^1. (18) 

-Tv/2 

The function !?((/') describes the physical properties of the tunneling barrier, and therefore depends on the details of 
the experimental configuration. For example, a D{(j)) = 6{4>) describes a "thick" barrier with dominant tunneling into 
the forward direction. Whereas a -0(0) = I/tt corresponds to a "random" barrier with equal probabilities of tunneling 
into all directions. 
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^max-bulk 

Fig. 8. The T = differential tunneling conductance as a function of bias voltage, eV, for a (110) surface and 
a (100) surface and three different sampling-cone sizes. An energy width (~ O.lAmax-buik) has been included to 
simulate experimental broadening. 

We present results for three types of model barriers. The barriers are described by a constant scattering probability 
into a "sampling cone" about the normal direction with an opening angle 6 cone- Specifically, we consider: (1) a 
thick barrier with Ocone — deg, (2) a random barrier with Ocone — 90 deg, and (3) an intermediate barrier with 
Ocone — 30 deg. In Fig. 8 we show the differential tunneling conductance for these three barrier models considering 
the two most natural surface to lattice orientations, i.e. a [100] surface {do = Odeg) and a [110] surface {(t>o — 45 deg). 
Note that for 0o = 45 deg one finds a considerable zero-bias anomaly [^,|| , while for 0o = deg one finds the bulk 
density of states averaged over the sampling cone. Thus one expects, for d-wave superconductors with ideal surfaces, 
a very different tunneling characteristic depending on surface orientation. We know of no experiment which has 
considered the tunneling spectrum as a function of surface orientation and which would therefore be in a position to 
either verify or refute these clear predictions of the d-wave model. 




0.0 0.2 0.4 0.6 0.8 1.0 

xl w 



Fig. 9. Typical order parameter profiles in the slab geometry when both a Bi and a B2 component are present. 
The data are for the (1,1,0) surface {(j>o — 45 deg) of the slab. The solid (dotted) lines correspond to the amplitudes 
of the Bi {B2) components. The curves are drawn for three different slab widths: lOvr^ (circles), Svr^ (squares), 
and S.Stt^ (diamonds). 
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B. Finite Slab 



We consider a slab of finite thickness, w, and infinite lateral extent, with two specular walls located at a; = and 
X = w. This is a very convenient setting for studying surface effects since a quasiparticle experiences the surface 
regions several times along its trajectory. This enhances the effects of surfaces on the density of states. This system 
was studied in detail by Nagato et. al. ||ll| for the case of a single-component i32-symmetry order parameter, and a 
surface to lattice tilt angle of (/)o = deg. They discussed the existence of a phase transition from the superconducting 
state to the normal state as a function of layer thickness. This scenario becomes more complicated, however, if we 
allow the mixing-in of sub-dominant symmetry components. In this case we find, for certain tilt angles (j)o, that a 
phase transition involving only one of the symmetry components can occur, leaving the system in a superconducting 
state of a different symmetry. In this way the normal state is not realized until one reaches much smaller sample sizes, 
if at all. 
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Fig. 10. The dependence of the order parameter components on the width of the slab: (a) The maximum of 
the Bi-symmetry component for various slab to lattice tilt angles, (po, and (b) the maximum of the _B2-symmetry 
component for the same tilt angles. Note that the system has a phase transition at a slab width w — 9.9^ for 
cpo = 45 deg. 

In order to exhibit these effects we consider the case of an isotropic cylindrical Fermi surface and compute the 
self-consistent order parameter in the same manner as discussed in [I]. We again assume that the Bi symmetry 
channel possesses the highest transition temperature, and allow for the mixing in of a B2 symmetry component 
with TcB2/Tc = 0.7. In Fig. 9 we plot a typical set of order parameter profiles for a surface to lattice orientation 
of (j)o = 45 deg and various slab widths. As the sample becomes thinner the Bi component is strongly suppressed, 
while the presence of the B2 component becomes enhanced. A plot of the maximum absolute value of each of the 
components as a function of the slab width is shown in Fig. 10 for several surface to lattice orientations. The angle 
(j}o = 45 deg is a special case since it represents the maximal pair-breaking orientation for the Bi component and the 
minimal pair-breaking (i.e. symmetric) orientation for the B2 component. Hence, the Bi component is favored in the 
bulk since it possesses the highest Tc, however the B2 component is favored at the wall for symmetry reasons. As 
a result, at a width w « 9.9^ we find a second order phase transition in which the Bi component abruptly vanishes 
leaving behind a spatially uniform i32-symmetry order parameter. As we move to angles away from 0o = 45 deg, we 
find that the phase transition is smeared out. 

In Fig. 11 we plot the superconducting DOS factor, in the vicinity of the phase transition for two different 
quasiparticle trajectories. The dotted lines show the form of v below the critical width and are just what one expects 
for a spatially uniform order parameter (i.e. they are identical to the bulk results). Above the critical width the 
spectrum exhibits a complicated band- like structure |l^]. Since we are considering a clean system with specular walls. 
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a quasiparticle will propagate without dissipation, reflecting an infinite number of times, always with the same angle 
of incidence. The local order parameter along such a trajectory has a periodic structure so that this scenario closely 
parallels the case of a particle in a periodic potential. The splitting of the excitation spectrum into bands is thus a 
consequence of Bloch's theorem. 



5.0 
4.0 
to 3.0 

H 

y 2.0 

1.0 

0.0 

4.0 

to 3.0 
H 

^ 2.0 
1.0 
0.0 



- x = 0^ 




(a) 


, J .■ 






J 


\ I 


(^=10° ('^^ 


w = 9.74§ 

w = 10.0^ 




„ / r = 0.7 " 







0.0 



0.5 

8/ A 



1.0 



1.5 



max-bulk 



Fig. 11. The superconducting density of states in a slab with a slab to lattice orientation (po ~ 45deg. Results 
are shown slightly above (it; — 10.0^) and below {w = 9.7^) the phase transition for two different trajectories: (a) 
(f) = Odeg, and (b) — lOdeg. 



We have also considered the case of a Si-symmetry order parameter coupled to a sub-dominant A2-symmetry 
component. In this case it is the tilt angle (j)o = 22.5 deg which is special, since it represents the symmetric orientation 
for the A2 component. Reminiscent of the previous case, we again observe a phase transition around w ~ 5.5^ in 
which the Bi component vanishes while the A2 component becomes spatially uniform. However, here the transition is 
first order and, moreover, several distinct near-by metastable states compete. In the transition region, the differences 
in free energy between the various states are so small that specific perturbations in any experimental configuration 
would surely be decisive. The superconducting DOS factor, exhibits distinguishing features on either side of the 
transition which, again, just parallels the previous case. 



IV. CONCLUSIONS 



We have studied the spectral properties of an anisotropically paired superconductor in the vicinity of a 
superconductor-insulator interface by considering two different physical geometries: a semi-infinite system, and a 
slab of finite width and infinite lateral extent. The stable bulk phase of the order parameter is always taken to trans- 
form like the Bi (dx'^-y^) representation of the D^h (tetragonal) group, however other attractive symmetry channels 
are also allowed to be present, which has a pronounced effect on the order parameter's behavior near the walls. The 
calculation of the self-consistent order parameter for systems such as these was discussed by the current authors in 
Ref. ]l3[ | and, we have carried over those results as input into the current calculations. All of our calculations were 
performed considering both an isotropic cylindrical Fermi surface, as well as a more realistic tight-binding Fermi 
surface; we find very little difference in the qualitative features of our results between the two Fermi surface models. 

The spectra at the surface of a semi-infinite system exhibit three characteristic features. First, one finds band-edge 
singularities at energies corresponding to the bulk gap magnitudes |Af,„;fc(p^. and |Af,„;A;(p^ q^^)! which yield 
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information about the gap anisotropy far from the wall. Secondly, one observes the formation of surface bound states 
which are of two distinct types. For any order parameter which changes its sign along the Fermi surface (regardless 
of its symmetry), one can always find situations in which a bound state will be formed at zero energy |qj20||. The 
formation of these states depends only on the relative sign of the asymptotic order parameter values (AhuikiPf [jj) 
and Af,uifc(P/ out)) along a quasiparticle trajectory. As a consequence, the mixing-in of sub-dominant symmetry 
components with the bulk (Bi) symmetry component has no impact on the presence of these states, since the order 
parameter is only altered within a few coherence lengths of the surface. Depending on the detailed nature of the 
interaction of the order parameter with the wall, however, a second species of bound states may appear at higher 
energies, if the order parameter is sufficiently suppressed in the surface region. The weight and location (in energy) 
of these bound states is directly related to the spatial dependence of the order parameter near the wall. Thus, one 
finds that the formation of these states can be correlated to both the strength and symmetry of the sub-dominant 
components. All of these features depend strongly on the surface to lattice orientation, and thus one can conceive of 
mapping out the full anisotropy of the order parameter by measuring the surface density of states for various surface 
orientations. 

The study of a d-wave order parameter in a slab geometry can be especially illuminating since the percentage of the 
sample volume which lies in the surface region is significantly enhanced. For the case of single d-wave order parameter 
in a finite slab, one finds phase transitions for certain surface to lattice orientations which results in a nontrivial 
temperature vs. width phase diagram |pl| . When the bulk Bi symmetry order parameter is allowed to couple to a 
B2 symmetry component, these phase transitions are washed out for most surface orientations. For certain "special" 
surface orientations however, phase transitions still occur (of either first or second order) which may take the sample 
into the normal state or even into a state with a different bulk symmetry. These special orientations are usually those 
in which the sub-dominant component does not pair-break. These phase transitions have a pronounced effect on the 
excitation spectrum since the order parameter always has a complicated spatial structure on one side of the transition 
while it is completely structureless on the other. 

The calculations presented in this paper were carried out for clean superconductors with perfectly reflecting surfaces. 
This is important to note because the presence of a diffuse boundary could result in significant qualitative effects. 
Indeed, our ongoing work indicates that for diffuse surfaces the distinction between the various surface to lattice 
orientations may no longer be as clear, since any one in-coming trajectory may mix with all out-going trajectories. 
This effect tends to smear-out, to variable degrees, some of the aforementioned phase transitions and characteristic 
features in the surface excitation spectrum. We conclude that a concerted effort should be made to fabricate "clean" 
surfaces, with specific surface to lattice orientations, if the maximal quantity of information on the anisotropy of the 
order parameter and the surface induced features in the density of states is to be extracted. 
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